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Abstract 

Time scales in turbulent and transitional flow provide a link between experimental 
data and modeling, both in terms of physical content and for quantitative assessment. 
The problem of interest here is the definition of time scales in an unsteady flow. Us- 
ing representative samples of data from GEAE low pressure turbine experiment in 
low speed research turbine facility with wake-induced transition, we document several 
methods to extract dominant frequencies, and compare the results. We show that 
conventional methods of time scale evaluation (based on autocorrelation functions and 
on Fourier spectra) and wavelet-based methods provide similar information when ap- 
plied to stationary signals. We also show the greater flexibility of the wavelet-based 
methods when dealing with intermittent or strongly modulated data, as are encoun- 
tered in transitioning boundary layers and in flows with unsteady forcing associated 
with wake passing. We define phase- averaged dominant frequencies that character- 
ize the turbulence associated with freestream conditions and with the passing wakes 
downstream of a rotor. The relevance of these results for modeling is discussed in the 
paper. 
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1 Introduction 


Dominant turbulent time scales characterize important physics of turbulence, and 
are featured prominently in models applicable to turbomachinery flows: the ratio k/s 
characterizes large eddies and is ubiquitous in the modeled expressions of the trans- 
port, production and destruction terms. The problem is that conventional methods 
(based on autocorrelation functions and Fourier analysis) are ill-suited to the treat- 
ment of unsteady data. In this paper, we compare the time scales obtained from 
conventional and wavelet tools, and illustrate unique capabilities of the wavelet algo- 
rithms downstream of a rotor. 

Three goals were addressed in the process. 

1. A quantitative comparison of integral scales and most-energetic, scales, based 
on conventional methods applicable to stationary signals (use of autocorrelation 
functions and Fourier spectra), shows that the two scales are not very different. 
This allows us to focus on most-energetic scales while allowing the results to be 
converted to the integral-scales favored by modelers. 

2. Then, we show that wavelet-based methods present an alternative to the Fourier 
spectra, with excellent quantitative agreement. Thus, the integral scale infor- 
mation can be obtained from wavelet-based results. 

3. Finally, we use the greater flexibility of the wavelet methods to extend the 
measurement of dominant time-scales to the unsteady environment of embedded 
stages. The flow characterization obtained in this manner may be compared 
with unsteady numerical results. 

1.1 Time scales and modeling 

Let us define modeling as the construction of systems simpler than the complex 
reality. The selection of relevant physics at the exclusion of unessential features is at 
the core of any scientific or engineering model, and sets the limits of its applicability. 
Indeed we might think of the compressible Navier-Stokes equations as ‘reality’, but 
the depiction of the fluid as a continuum is certainly a model, with known limits of 
accuracy. Thus, it is proper to envision a hierarchy of models of reality, with high- 
level models providing a point of reference in terms of accessibility and accuracy of 
the solutions. Lower level models are obtained by ignoring some terms, or by using 
simpler alternatives for analytically-difhcult terms; explicit reference to the higher- 
level model may even be eliminated in some cases. 

An engineering model will be conceived and evaluated based on considerations of 
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• physical content: qualitative or quantitative model predictions can be expected 
only if the relevant dynamics are included in some form. It would make little 
sense to ignore turbulent transport in studying flow through turbomachinery, 
whereas the same assumption is acceptable in the context of homogeneous tur- 
bulence. 

• analytical form: with no known general form of the transport terms in averaged 
equations for turbulent flows, eddy- viscosity models are frequently used. 

• quantitative parameter tuning: by adjusting model constants, model predictions 
can have quantitative value to the engineer. 

The material presented in this report might affect physical content as well as quanti- 
tative adjustments in turbomachinery flow modeling. 

The selection of physical content in a model is usually based, more or less explicitly, 
on scale considerations. Molecular and sub-molecular dynamics are not considered 
essential is describing flows over scales much larger than the mean-free-path; climatic 
drift is not essential in updating weather predictions for the next few days; and so 
forth. In some circumstances, the energy spectrum characteristic of the phenomenon 
will exhibit gaps, which effectively isolate some spectral ranges and allow us to study 
them in isolation; in other cases such as turbulence, broad-band spectral activity 
dictates a more cautious approach to model construction. 

In this introduction, we will use length- or time-scales indifferently; the data 
analysis will focus on the more directly accessible time scales. Based on the features 
of a typical power spectrum of a turbulent flow as on Fig.l 1 , one can identify several 
scales: starting from low-frequencies, 

1. the point, at approximately 1 kHz, where the energy content starts decreasing 
is distinctive (macroscale), 

2. as is the range of scales, around 4 kHz, in which power-law behavior is exhibited, 

3. as is the point, approximately 8 kHz, at which the spectrum tails down (mi- 
croscale or limit of resolution). 

Appropriate definitions help us pinpoint these scales: integral scales, Taylor and 
Kolmogorov microscales, etc. The interest of these scales in modeling is that they are 
representative of certain dynamical effects at the exclusion of others, in spite of the 
absence of a spectral gap. For instance, 

1 it is immaterial for this illustration that the high frequency drop-off is actually due to low-pass 
filtering rather than to actual viscous dissipation 
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X— wire: rxw324a02 - inlet to N1 

u and v power spectra 

Figure 1: Longitudinal (dotted line) and transverse (solid line) power spectral density 
at inlet of nozzle 1. 

1. the Kolmogorov microscale is characteristic of the viscous dissipation of energy, 
and (at large Reynolds numbers) contains very little energy and accomplishes 
very little transport or mixing. 

2. the Taylor microscale is representative of the absence of either viscous or trans- 
port effects: from the equations of motion, only unsteady inertia effects (con- 
vection/pressure/spectral transfer) remain. 

3. the macroscales can be defined as either integral scales or as energy- containing 
scales: they are representative of large eddies that accomplish much of the 
mixing in the turbulent flow, and of the production of turbulence energy, in 
isolation from significant viscous effects. The remainder of this paper will focus 
on the macroscales. 

This scale separation of the various dynamical terms is at the basis of much turbulence 
modeling. For example, Kolmogorov’s equilibrium range model (small scale) combines 
inertial and viscous effects. Closer to the subject of this report, in modeling the 
transport terms in a turbulent flow, macroscale parameters should be used. Using 
viscous scaling for turbulent transport would generally be regarded as bad modeling 
on account of the inappropriate physical content. Thus, a suitable understanding of 
the turbulent macroscales and their quantitative evaluation are at the core of any 
model of turbulent transport. 
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Time scales are particularly indicated in the context of this project because of the 
following reasons: 

• the data was collected as time series, so that characteristic frequencies are more 
natural than a corresponding wavenumber - the conversion between the two 
makes use of Taylor’s Hypothesis. 

• it is generally accepted (see e.g. Tennekes & Lumley, p.70) that the ratio of 
time scales between two phenomena is a measure of how well correlated they 
are and of how effectively they exchange energy. For example, small eddies 
behave independently of the large scale turbulence; a process is adiabatic if it is 
fast relative to applicable time scales of energy transfer; etc. Thus, quantitative 
evaluations of time scales will provide a measure of physical interactions within 
the flow. 

• most turbulence models (e.g. 2-equation models) provide ready access to large- 
eddy time-scales, in the form of k/s found in the construction of transport, 
production and destruction terms. 

1.2 The data 

In this study, we focus on the time-scales of anemometry data, in an attempt to 
identify sensitive algorithms capable of characterizing the physics of transitioning 
flows subject to rapid straining and unsteady freestream turbulence. The times scales 
are one channel by which experimental information is normally reflected in model 
results, and we expand the applicability of the idea to unsteady flows. The data 
is from an experiment conducted by Dr. David Halstead at GE Aircraft Engines 
(GEAE). The experiment was conducted at GEAE Low Speed Research Turbine rig 
on a two-stage low-pressure turbine. The configuration is shown in Fig. 2. Halstead 
took measurements with X hot-film probe at the inlet and exit of each blade row, 
with hot-wire probe in the boundary layers of the first and second stage nozzles, 
and with hot-film arrays surface-mounted on the second stage nozzle and rotor. The 
complete data set is the object of extensive analysis presented by Halstead et ah 
[4, 5]. The work presented in this report is limited to representative X-probe data 
taken at the inlet and the exit of the first stage Nozzle and Rotor, designated as N1 
and Rl, respectively. Work utilizing quasi-shear stress data from ting surface hot-film 
arrays was performedby Lewalle [12]. It is expected that the methodology and data 
processing presented here will lead to further method development and extensive data 
analysis that will advance modeling of unsteady flows in turbomachinery. 
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Figure 2: Schematic of the facility 

1.2.1 Inlet conditions 

At the inlet, unspecified ambient turbulence is present. On Fig. 3, we show an 
instantaneous trace, an ensemble average, the fluctuation around that average, and 
the skew at the N1 inlet. There is no appreciable periodic component of the incoming 
turbulence. As a result, the ensemble mean velocity reflects mostly the random 
fluctuations of the incoming flow, with an mis value presumably decreasing as the 
square root of the number of traces in the ensemble (Mean Value Theorem). Indeed 
we see that the mis of the ensemble mean is smaller than the mis of each trace by 
one order of magnitude. We will study the fluctuations around this mean. 


1.2.2 Rotor exit conditions 

On Fig. 4, we show an instantaneous trace, an ensemble average, the fluctuation 
around that average, and the skew. In contrast to the N1 inlet data, a strong periodic 
component can be identified on the raw data, and the ensemble mean shows strong, 
smooth periodicity with mis comparable to that of the raw signals. Minor variations 
from period to period in the ensemble mean reflect tint turbulent fluctuations, which 
have a similar mis. 

The ensemble mean velocity components were subtracted from the processed 
traces, using the formulae listed in Halstead [5]. This enables the ‘random’ tur- 
bulence to stand out from the periodicity of the wake shedding from Rl. Then we 
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Figure 4: Streamwise velocity at the R1 exit. From bottom, to top: raw data, ensemble 
average (300 traces), fluctuation around the average, and .skew. 
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reproduced the statistics provided by Halstead, as a test of integrity of the processed 
d at a. 

1.3 Overview 

We proceeded with the calculation of several time scales from this data. 

• calculation of the integral time scales, obtained by integration of the autocor- 
relation functions for the streamwise and transverse velocity components. 

• calculation of the integral time scale from the zero-frequency estimate of the 
1-dimensional power spectrum 

• identification of the most-energetic scales from the Fourier power spectral den- 
sity. 

• identification of the most-active transporting scales from the Fourier cross- 
spectral density. 

• identification of the most-energetic scales from the Morlet wavelet mean power 
spectrum. 

• identification of the most-energetic scales from the Mexhat wavelet mean power 
spectrum 

• identification of the most-active scales from the Mexhat Reynolds stress decom- 
position 

Detailed presentation of these methods can be found in Hinze [7], Lewalle and Ash- 
pis [11], and other references itemized in the Appendices. We compare and interpret 
the results as follows: 

• at the inlet section, we tabulate the various statistics and note differences. We 
compare our results to classical isotropic-turbulence relations between integral 
scales and energy spectra [7]. This establishes the validity of the statistics for 
parameterization and characterization purposes. 

• at the embedded station, we confirm that mean properties can be calculated 
by any of the same methods, and note that the mean statistics are .father 
unsatisfactory as descriptors of the complexity of the flow. 

• the flexibility of wavelet spectra in relation with conditional sampling is illus- 
trated with the construction of phase-dependent time scales. The variations in 
dominant turbulent scales during one period of wakq-passing are mapped and 
discussed. 
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2 Method Validation on Stationary Data 

We first document the time scales for stationary turbulence at the entrance to the 
first rotor (Nl). The purpose is to show that similar results are obtained from integral 
time scales, dominant scales of the Fourier spectrum, and dominant wavelet scales. 
Unsteady data at the rotor exit (Rl) will be treated in the next Section. 

For brevity of the presentation, only essential formulae will be listed in this section. 
Details of implementation, differences between continuous and discrete formulations, 
and specifics relating to factors of 2 tt in the formulae, are listed in Appendix 1 for 
reference. 

2.1 Integral time scales 

Let us focus first on one velocity signal, say u(t). Its standard deviation is measured 
by u'\ 

u' 2 = lim ^ / \ u{t)\ 2 dt. (1) 

co Zl J-T 

The longitudinal autocorrelation function for a stationary signal u(t) is defined as 

u' 2 f u {t) = u{T)u{t + r) (2) 

The autocorrelation functions measure the gradual loss of correlation as the time lag 
between measurements of the given component u increases. Analytical and exper- 
imental properties of these functions can be found in standard references [7]. For 
diagnostic purposes, it is convenient to summarize the information about the drop-off 
by a single number: the longitudinal integral time scale is defined as 

p oo 

T u ,int = / fu{t) dt , (3) 

Jo 

In addition to the longitudinal autocorrelation function (Ecp 2) and integral time 
scale (Ecp 3), we define the transverse autocorrelation function 

v / 2 g v {t) = v{r)v{t + t) (4) 


and transverse integral scale 

poo 

T v , int = g v {t) dt. (5) 

Jo 

Both autocorrelations are shown on Fig. 5. The correspondence between experimental 
results and model parameters is normally based on these integral scales. Integral 
scales are clearly defined, can be computed rapidly, and capture the ‘size* of the large 
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(energy- containing) eddies. The natural time scale occurring in relation to these large 
eddies (a measure of eddy turnover time) in k-s models is the integral time scale 

Tu,int = - (6) 

£ 

(see e.g. Wilcox, p.90-92), and all model constants are adjusted accordingly. The 
values can be converted to integral length scales using Taylor’s Hypothesis under a 
convection velocity U: 

A / = UT Ut int (7) 

To make comparisons possible with results from later Sections, we had to choose 
between time scales or frequencies, and settled on the latter. The frequencies corre- 
sponding to the time scales above were calculated as 

fu, int = = 0.606A’i7~ (8) 

1 u,int 

(respectively) 

fv, int = = 1.106 kHz. (9) 

1 v,int 

We note that the approximate factor 2 between these two scales agrees with the 
theoretical relation between the longitudinal and transverse integral scales (Hinze [7], 
Eqs. 3-74 and 3-75) for isotropic turbulences 

Practically, the computation of the autocorrelation function is very time con- 
suming. Furthermore the slow convergence of the function for large times makes its 
numerical integration rather inaccurate. For these reasons, an alternative definition 
based on the Fourier spectrum Ep (see Fig. 1) is generally preferred. It is easy to 
show (e.g. Hinze [7], p65) that the integral scale defined above can be calculated from 
the formula 

T _ Ef( 0) 

4 u ,z 

The convergence problem for the autocorrelation function now takes the form of 
inaccuracies in the calculation of the low-frequency part of the spectrum. However, as 
the energy density is constant in this region of the spectrum, extrapolation of any well- 
converged part of the spectrum to lower frequencies gives reliable estimates of Ep( 0). 
Thanks to the use of FFT’s, the integral scales can be calculated very efficiently. 
These observations combine to make this the preferred method of calculation of an 
integral scale. The results are as follows: 



fu, Eo = o.mkHz (ii) 

and 

f v , Eo = 1-387 kHz (12) 
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X-wire: rxwu24a02 — 'nlet to Ml 
u — a utocorrelation function 



X — wire: rxw324a02 — inlet to N 1 
v— autocorrelation function 


Figure 5: Longitudinal (top) and transverse (bottom) autocorrelation at the inlet of 
nozzle 1. 
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ill good agreement with the values listed above. The 10% to 20%) discrepancy is 
accounted for by incomplete convergence in the tail of the autocorrelation function 
and in the low-frequency portion of the spectrum. 

2.2 Dominant Frequencies 

An alternative characteristic time scale based on the Fourier spectrum is the domi- 
nant (most energetic) spectral range. Formulae are reproduced in Appendix 2. It is 
generally known that the most energetic frequency is identified as the frequency cor- 
responding to the maximum of the f Ep(f) curve. Computations of this statistic are 
very fast, and accuracy is good as long as the power spectrum shows little scatter. If 
necessary, smoothing of the power spectrum can be introduced to mitigate the effect 
of spurious spikes. 

Although the definitions of the integral scale (and its corresponding frequency) 
and of the dominant frequency are not identical, they are related in the sense that the 
turbulent eddies most closely associated with the integral time scale are the ‘large’ 
eddies - the ‘ energy-containing ’ eddies. Thus, the integral scalq and the dominant 
scale characterize the same spectral range, and the numerical values should be of the 
same order of magnitude. Since model constants are adjusted based on the values 
of integral scales, a quantitative conversion factor from dominant to integral scales 
is required if the dominant scale is to be used. In appendix 3, a simple analytical 
model of the power spectrum provides a possible conversion factor. This issue will be 
addressed further in the conclusions. 

The Fourier power spectra for u and v are shown on Fig. 1. The spectra are 
generally featureless, but the two components are appreciably different. In contrast 
to the integral scales, a relevant time scale can be defined as a local (pointwise) 

property of the spectrum, i.e. the scale corresponding to the most energetic eddies. 

The dominant frequency defined in this manner is easily identified at the peak of 
f.E(f), as shown on Fig. 6 

The experimental scatter in the spectra will lead to some minor differences in the 
determination of the dominant frequency. From Fig. 6, we get the values 

/eu = 0.666 kHz (13) 

and 

f Ev = 1.188kHz. (14) 

The values are in close agreement with the integral scales provided above. 
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Figure 6: Longitudinal (dotted line) and transverse (solid line) energy content at inlet 
of nozzle 1. 

2.3 Wavelet-based dominant scales 

Wavelets have emerged in recent years as a useful alternative to Fourier analysis for 
the case of unsteady data. Some fundamental relations about wavelet analysis are 
listed in Appendix 4. We apply the wavelet processing to the stationary data to 
demonstrate the equivalence with Fourier methods in this case. Given this reliability 
check, the justification for using wavelet methods will be found in the analysis of 
unsteady data in the next Section. 

We used two different wavelets in this analysis (see Appendix 4 for details). The 
Morlet wavelet is a complex- valued wavelet, most appropriate when local periodicity 
is expected. The Mexican Hat wavelet is real, and is better suited to identify isolated 
features in the signal. Temporal localization is best viewed from the time series; spec- 
tral features are clearly visible on the Fourier transform. Between these extremes, a 
combination of temporal and spectral information is captured by wavelets, and differ- 
ent wavelets correspond to a different emphasis on temporal vs. spectral resolution. 
In this framework, the Morlet wavelet achieves better spectral resolution with poor 
localization, whereas the Mexican Hat wavelet favors temporal resolution at the cost 
of some spectral resolution. This point is illustrated by the respective power spectra, 
in which the smoothness of the Morlet spectrum is intermediate between the Fourier 
and Mexican Hat spectra (see Figs. 7 and 8). The comparison with Fig. 6 shows 
excellent agreement between the three methods. 

The most energetic frequencies are identified from these figures as 

fu-Mor = 0.7207-HA, (15) 



1000 1QQ0O 

Frequency f (Hz) 


X — wire: rxw324a02 — inlet to N1 
u and v Fourier energy 
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X-wire: rxw324a02 — inlet to N1 
u and v Morlet energy 


Figure 7: Streamwise (dotted line) and transverse (solid line) Morlet energy content 
at inlet of nozzle 1. 



X-wire: rxw324a02 — inlet to N1 
u and v Morlet energy 

Figure 8: Stream wise (dotted line) and transverse (solid line) Mexhat energy content 
at inlet of nozzle 1. 
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X-wire: rxw324a02 — inlet to N1 

.... Mexhat and Fourier spectral stress 

Figure 9: Fourier- (solid line) and Mexhat (dashed line) spectral stress distribution. 


fv — M or = 1.079 kFdz, 

(16) 

fu-Mex = 0.796 kHz 

(17) 

fv-Mex = 1.321 kHz. 

(18) 


We conclude that the time scales calculated by conventional methods and wavelet 
results are in agreement. 


2.4 Spectral distribution of the turbulent stress 

X-probe data enables one to calculate scales relevant to the turbulent transport. Mul- 
tiplication of the co-spectrum by the frequency yields the spectral stress distribution, 
shown on Fig. 9 

A similar decomposition of the Reynolds stress was presented in AIAA-94-2323, 
and further analyzed in NASA-TP-3555. The Mexhat stress distribution is also shown 
on Fig. 9. We note again the general agreement between the two methods. The 
dominant frequencies are identified, respectively as 

fcoFou = 0.924kHz (19) 

and 

f co M ex = 0.975 kHz. <21)1 
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2.5 Comparison and interpretation 

To summarize, the numerical values of characteristic or dominant frequencies calcu- 
lated in the previous subsections are reproduced in tabular form. 

Table 1: Comparative table of dominant frequencies at the N1 inlet. 


fu, int 

fv, int 

fu, Eo 

fv, Eo 

0.606 kHz 
1.106 kHz 
0.671 kHz 
1.387 kHz 

f Eu 

0.660 kHz 

fEv 

1.188 kHz 

f coFou 

0.924 kHz 

fu—Mor 

0.720 kHz 

fv—Mor 

1.079 kHz 

f u—Mex 

0.796 kHz 

fv—Mex 

1.321 kHz 

fcoMex 

0.975 kHz 


The numerical values of possible relevant frequencies fall into three main cate- 
gories. 

1. Longitudinal (strearnwise) scales fall in the range of 0.606A’JL~ to 0.796 kHz. The 
lower value suffers the most from known problems associated with the accuracy 
in the tail of the autocorrelation function, and may be ignored. The frequency 
0.671 kHz should be preferred in relation to the integral scale. Energy- dominant 
scales are then 0.660A’i7~ for the Fourier spectrum (affected by the spurious 
peaks resulting from incomplete convergence), 0.720A’i7~ from the. Morlet spec- 
trum and 0.796 kHz form the Mexican Hat spectrum. 

2. Transverse components yield larger values. We already noted that, in isotropic 
turbulence, a factor of 2 would be expected, so the inlet turbulenc# results 
are consistent with nearly-isotropic. turbulence. The most reliable value from 
autocorrelation is 1.387 kHz, while the Mexican Hat dominant frequency is 
1.321kHz, a 5% discrepancy. 

3. Finally, values of the frequency contributing the largest part of the Reynolds 
stress is intermediate between the above. 

We conclude that 
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X-wire:, ■|®xw-324e04 — inlet to N2 
u — autocorrelation i i notion 


Figure 10: Longitudinal autocorrelation at the inlet of the second stage (N2). 

1. while the Reynolds stress scales are consistent with the u- and v-scales, their 
use is subject to data availability and presents no obvious advantage over the 
1-component alternatives. 

2. differences between the u- and v-components are expected from the theory of 
isotropic turbulence. 

3. the integral scale method and wavelet method yield results within 7% of each 
other. 

Having shown the equivalence of the methods when applied to stationary turbu- 
lence, we now turn to the unsteady case. 


3 Unsteady Flow at the R1 Exit Section 

3.1 Further Validation for Mean Time Scales 

The streamwise (longitudinal) autocorrelation function is shown on Fig. 10. We 
note the negative values of the autocorrelation within time lags of approximately .3 
ms. This overshoot leads to the partial cancelation of the integral (Ecp 3). This 
problem is combined with statistical fluctuations in the ‘tail’ (large time lags) of the 
autocorrelation, and leads to inaccurate values of l lie integral scale. This problem 
is alleviated by the extrapolation of E\ to low frequencies; only these values will bq 
tabulated below. 


NASA/TM— 2004-209452 


18 




The Fourier energy distributions for streamwise and transverse fluctuations are 
shown on Fig. 11. The N1 results (Fig. 6) are shown to scale. The relative magnitudes 
of the incoming and wake-induced turbulence can be observed, as is the shift toward 
higher frequencies. However, the separation in scale is not sufficient to identify the 
scales as distinct on thes# figures. We see a general shift toward higher frequencies, 
particularly in the case of the transverse component. This result is true also for the 
Morlet- and Mexhat energy distributions, not shown here. 

A different picture emerges from the spectral distributions of Reynolds stress. Fig. 
12 shows the Fourier- and Mexhat versions. Beside the good agreement between the 
two methods, the changes from the N1 inlet are noteworthy. The stress distributions 
are in fact double-peaked, with a low-frequency peak identified at 500 Hz (the wake- 
passing frequency, in spite of the subtraction of the ensemble mean), separated from 
a weaker peak in the 5 kHz range by a spectral range of weakly negative correlations. 
Further study will determine how spectral regions can be associated more ; clearly with 
the wake passings. 

As for the N1 data, the characteristic frequencies are tabulated below for easier 
comparison. 

Table 2: Comparative table of dominant frequencies at the N2 inlet. 


o o 

1.254 kHz 
0.991kHz 

f Eu 

1.122 kHz 

f Ev 

1.023 kHz 

f coFou 

.528 kHz 

fu—Mor 

1.079 kHz 

fv—Mor 

.975 kHz 

f u—Mex 

1.390 kHz 

fv—Mex 

1.079 kHz 

f coMex 

.559 kHz 


3.2 Comparison and interpretation 

The first observation on Table 2 is that the 2:1 ratio of longitudinal to transverse 
scales predicted by l lie theory of isotropic turbulence no longer holds true v Second, 
the methods based on single components of velocity agree with each other. Finally, 
the dominant scale in the spectral distribution of Reynolds stresses is surprisingly 
low, and very close to the wake passing frequency. 
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X-wire: rxw324c04 - - inlet to N2 
u— Fourier power spectrum 



X-wire: rxw324c04 - - inlet to N2 
v— Fourier power spectrum 


Figure 11: Streamwise (top) and transverse (bottom) energy content at inlet of nozzle 
2. The nozzle I results are reproduced in dotted line. 
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100 


10000 


1000 

Frequency (Hz) 


X — wire: rxw324c04 — inlet to N2 

.... Mexhat and Fourier spectral stress 

Figure 12: Spectral stress distribution: Fourier (solid line) and Mexhat wavelet (dotted 
line). 
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Even if we overlooked this last point, the characterization of the data is as yet 
unsatisfactory. The spectral stress distribution can be taken as an example (Fig. 
12). Three scales are identified — two associated with positive correlations of u and 
v, separated by a negatively correlated patch of scales. It is clear that several flow 
phenomena are at play, but the time averages do not allow us to separate them. Given 
the flow conditions, we surmised that the different spectral ranges are associated with 
the unsteady wake-passing. A phase conditioning similar to that reported in Schobeiri 
et al. [14] is indicated, and is the object of the following section. 

3.3 Conditional Scales from Wavelet Processing 

One of the dominant features in the preceding sections is the overall agreement be- 
tween the various (Fourier vs wavelet) methods. However, it is clear by inspection of 
the raw data, and also from the wavelet maps shown in Section 2, that the character- 
istics of the fluctuations are modulated by the wake passings. To associate a unique 
frequency with these fluctuations is bound to mask some of the flow physics. 

The key to the determination of local time scales is found in the wavelet methods. 
A variant of the conditional sampling method introduced in Higuchi et al. [8] and 
used in Schobeiri et al. [14], was developed for this study. It consists of calculating 
the norm-Morlet transform 2 of the ensemble mean, hireling its dominant frequency 
/ + , and calculating its phase </>(/ + ,t) at that frequency and as a function of time; 
see Walker et al. [15] for details and examples. Then, the time-integral of energy in 
the Parseval relation (Ecp 24) is restricted to those contributions that coincide with 
a narrow range of phases: essentially, we construct phase-locked energy distributions, 
i.e. the signal’s energy content as a function of frequency and phase relative to the 
ensemble mean. 

The algorithm is as follows. The ensemble average velocity was constructed, and 
its dominant frequency (500 Hz, corresponding to the blade passing frequency) de- 
termined. Then, the phase-Morlet transform gave the phase of the ensemble average 
as a function of tim@y For a sinusoidal signal, the phase history would consist of a 
straight line increasing by 2 tt per period, with periodic resets (mod(27r)). This yields 
a uniform probability distribution of the phase angle of the interval (0, 2? r). The phase 
distribution for the ensemble velocity (streamwise, say) at the R1 exit is very nearly 
uniform. In the expressions of the wavelet energy and stress distributions (Eqs. 45 

2 From Eqs. 42 and 43, it is evident that the Morlet, wavelet transform is a complex function of 
the real variables k and t. The norm of this function contains the energy information (see Parseval’s 
theorem, above), while local phase information is defined by the polar angle between the real and 
imaginary parts: 

u g (f,t) = u norm (f,t)e' 4if ’ t] (21) 
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100 1000 10000 

Frequency f (Hz) Ens. mean vel 


X-wire: rxw324a02 — inlet to N1 
u — Mexhat power spectrum 

Figure 13: Streamwise frequency/phase Mexhat energy distributions at the N1 inlet. 
Morlet statistics are very similar. 

and 48), the time integral can be restricted to a given phase, and the distributions 
plotted as a function of this phase. The result is examined in this section. 

At the N1 inlet section, where there are no wakes and the ensemble average velocity 
is extremely small, the dominant scales can be expected to be independent of the 
phase. Indeed, Fig. 13 shows the frequency /phase energy content of the signal, 
with the dominant time scale at approximately 700 Hz corresponding to the averages. 
Similar results are obtained for the transverse component, using either Morlet or 
Mexhat wavelets; and also for the the Reynolds stress distribution. 

These results are in sharp contrast to the Rl exit data. Here, the wakes of N1 
(chopped by the Rl) and Rl pass at the measuring station periodically, and are 
associated with very different dominant scales. 

Fig 14 shows the frequency /phase energy content of the fluctuations. For the 
streamwise component, a roughly uniform phase distribution corresponding to the 
average peak is observed around 1 kHz. More energetic fluctuations at a higher 
frequency (3 to 5 kHz) are observed around <f> = 2. By comparison with the ensemble 
mean signals (Fig. 7), we associate this smaller-scale turbulence with the rotor wakes. 
The transverse component shows stronger phase dependence of the energy content. 
The wake turbulence is now seen at higher frequencies (5 kHz), while the dominant 
mean peak is seen to result from contributions from phase angles larger than it. 

A slightly different picture emerges from the Mexhat analysis. Fig. 15 shows a 
very strong phase localization of the streamwise Energy distribution around <f> = 2, 
with weaker fluctuations for the remainder of the period. In contrast, the transverse 
component shows a nearly uniform (mean) distribution for all phases, and an isolated 
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Frequency f (Hz) Ens. mean vel 


X-wire: rxw324c04 — inlet to N2 
u Morlet power spectrum 



Frequency f (Hz) Ens. mean vel 

X-wire: rxw324c04 — inlet to N2 
v— Morlet power spectrum 

Figure 14: Streamwise (top) and transverse (bottom.) frequency/phase Morlet energy 
distributions at the R1 exit. 
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X-wire: rxw324c04 — inlet to N2 
u — Mexhat power spectrum 



Phase relative to ensemble mean 

X-wire: rxw324c04 — inlet to N2 
v— Mexhat power spectrum 

Figure 16: Streamwise (top) and transverse (bottom) most-energetic-frequency as a 
f unction of phase at the R1 exit. 
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X — vv ire: rxw324cQ4 — inlet to 112 
Mexhat R-stress spectrum 

Figure 17: Frequency/phase Mexhat stress distributions at the R1 exit. 

island at much higher frequency at the phase of the wake passing. 

The most energetic frequency (the ‘ridge’ on Fig. 15) is plotted on Fig. 16. We 
see a doubling of the dominant frequency at <f> = 2, for both components of velocity. 
The peak is more sharply defined for the transverse component. The streamwise 
component shows another increase (missing from the v-data) out-of phase with the 
wake-passing by tt/2. 

Finally, the stress distribution is shown on Fig. 17. We see the sign reversal of 
the stress almost exactly at <f> = 2. We interpret this reversal as matching the change 
in sign of the ensemble mean velocity gradients across the wake. We also see that the 
dominant scales on either side of the wake are not equal. The larger scales at <j> > 2 
could be associated with the , N1 wakes, and the yet larger ones at <j> < 1 remain to 
I >e interpreted. 

The next series of figures (18 and 19) show the phase-dependent energy content, 
and the dominant frequencies, for the raw data, including the ensemble velocities. 
Although the dominant frequencies are slightly lower than with the ensemble velocity 
subtracted, the wake passing frequency (about 500 Hz) is not dominant. 

3.4 Discussion 

The calculation of the frequency /phase energy and stress plots relies on the condi- 
tional sampling of the energy maps, i.e. the time-frequency distribution of energy. It 
may be possible to achieve this result with windowed Fourier transforms, but wavelet 
transforms are the natural technique to reach this goal. Having ascertained that 
Fourier and wavelet methods given essentially the same results with the more sta- 
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X-wire: rxw324c04 — inlet to N2 
Mexhat R — stress spectrum 


Figure 19: Frequency /phase Mexhat stress distributions at the R1 exit. 

tionary signals at the N1 inlet, we are justified in using the better-suited wavelet tools 
at the N2 inlet. 

Morlet and Mexhat results, both for energy and Reynolds stress, agree that the 
dominant scales vary with the phase of the wake-dominated ensemble mean velocity, 
and that the scales are appreciably shorter at a phase <j> ~ 2 corresponding to the 
smallest values of the streamwise ensemble velocity, i.e. the wake passings. There- 
fore, we conclude that the frequency /phase plots discriminate quantitatively between 
freestream and wake turbulence, their energy content and dominant frequencies — 
parameters at the core of k-s unsteady models. 

However, there are differences in the particulars of the distribution of energy 
in the frequency /phase plane depending on the wavelet. This is a reminder that 
spectral methods project the signal on a template: infinite complex sin# wav# in 
Fourier transforms, a wave packet in Morlet wavelet transforms, a single bump in 
the case of Mexhat wavelet transforms. The key underlying issue is the trade-off 
between frequency and time resolution — heuristically associated with Heisenberg’s 
uncertainty principle. The progression from Fourier to Morlet (with decreasing values 
of ~ 0 ) to Mexhat corresponds to increasingly accurate temporal localization of ‘events’, 
at the expense of a deteriorating spectral accuracy. The user should choos# the best 
tool accordingly. In our case, we see no reason to assume even local periodicity of the 
fluctuations around the ensemble mean, and the Morlet template appears somewhat 
contrived (an unrealistic model for local fluctuations); we view the Mexhat results as 
most informative. 
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4 Conclusions 

In this report, we have accomplished three things: 

1. We have developed flexible alternatives to the integral time scales as a way 
to characterize unsteady data. We have shown that the various definitions 
yield comparable results on modeled spectra. We also show agreement between 
integral scales and dominant scales for all data sets available for this preliminary 
study. 

2. Having established that wavelet tools yield results equivalent to Fourier meth- 
ods, we used the time-frequency resolution of wavelets to map phase-dependent 
properties. We were able to quantify spectral energy content and dominant 
frequencies as a function of plia>e. Thus, the effect of wake passing at the N2 
inlet could be characterized. 

The experimental data made available to us by GE are very rich in flow physics 
pertaining to wak|miduced transition on compressor blades. An obvious application 
of this information in the construction and testing of improved models for unsteady 
transition, and the comparison of model and experimental results. 

Given the reliance of modelers on the integral scale, it would be desirable to pro- 
vide a conversion factor between the integral and the dominant scales. As expected, 
the ratio of these scales turns out to be of the order of 1 in all cases where both can 
be calculated. A more detailed comparison leads to the following table: 


Table 3: Conversion factor from integral to dominant frequencies. These values should 
be compared to the ~ .75 number obtained in Appendix 3. 


N1 fu—Mor/fu,Eo 

N1 fv—Mor/fu,Eo 

fu—Mor/fu,Eo 
fv—Mor/fu,Eo 

0.720 / 0.671 = 1.07 
1.079 / 1.387 = 0.78 
1.079 / 1.254 = 0.86 
0.975 / 0.991 = 0.98 

N1 fu—Mex/ fu,Eo 

N1 fv—Mexf fu,Eo 

fu—Mex / fu,Eo 
fv—Mexf fu,Eo 

0.796 / 0.671 = 1.18 

1.321 / 1.387 = 0.95 

1.390 / 1.254 = 1.11 

1.079 / 0.991 = 1.09 


We note that the Mexhat statistics tend to be slightly larger than comparable 
Morlet statistics, and that all values are larger than the .75 expected from the Pao 
model in Appendix 3. It is apparent that the particulars of the spectral distribution 
of energy, and the particulars of the analyzing wavelet, introduce some variability 
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ill the quantitative result. This variability is inherent in the statistical reduction of 
complex information to simple numbers. The real question lies in the impact this 
variability may have on modeling. Experience will tell what adjustments to model 
constants, if any, are required to ensure agreement between wavelet-dominant scales 
and modeling outcomes. The fact remains that, in the unsteady environments en- 
countered in turbomachinery and other complex flows, conventional definitions cease 
to be meaningful; this report is a first step in devising alternatives. 
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Appendix 1 — Integral Scales 


Let us focus first on one velocity signal, say u(t). Its standard deviation is measured 
by u' (see Eq. 1). One of the relevant time scales is characteristic of the persistence 
of correlations between lagged signals, measured by the autocorrelation function. 
The longitudinal autocorrelation function for a stationary signal u(t) is defined by 
Ecp 2. The autocorrelation functions measure the gradual loss of correlation as the 
time lag between measurements of the given component u increases. Analytical and 
experimental properties of these functions can be found in standard references [7]. 
For diagnostic purposes, it is convenient to summarize the information about the 
drop-off by a single number: the longitudinal integral time scale is defined by Ecp 3. 
Practically, the computation of the autocorrelation function is very time consuming. 
Furthermore the slow convergence of the function for large times makes its numerical 
integration rather inaccurate. For these reasons, an alternative definition is generally 
preferred. Start from the Fourier transform u{f) of the signal u(t): 

/ OO 

u{t)e-' 2iKft dt (22) 

-CO 

Then, the Parseval theorem, 

/ CO poo 

\ u(t)\ 2 dt= df | ii(f) | 2 , (23) 

-co J — CO 

leads to the definition of the power spectrum. Following Hinze (p.57), we define 
Ep{f) by the relation 

p CO 

u 2 = E F {f)df, (24) 

Jo 

Then it is easy to show (e.g. Hinze [7], p65) that the integral scale defined above can 
be calculated from Ecp 10. The convergence problem for the autocorrelation function 
now takes the form of inaccuracies in the calculation of the low-frequency part of the 
spectrum. However, as the energy density is constant in this region of the spectrum, 
extrapolation of any well-converged part region of the spectrum, extrapolation of 
any well- converged part of the spectrum to lower frequencies gives reliable estimates 
of Ep{ 0). Thanks to the use of FFT’s, the integral scales can b#: calculated very 
efficiently. 

In the case of digital time series, the formulae are adapted as follows. We denote 
by u i , i = 0, ...(N — 1) the N values of velocity collected at time intervals A so that 
ti = (i — 1)A. It is assumed that ip has zero mean. Then 


u 


12 


1 

N 


N - 1 


£ 


(25) 
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The autocorrelation function is defined as 


J i 


1 M—l 

j=o 


(26) 


where M and the range of f s are restricted to match the number of data points. For 
example, if M = N/2, f t will ’ be defined over the range i = 1, (N / 2 — 1). Then, 
the integral scale is given by 


T ■ — 

- J - u,tnt 


2A 

w 


N/ 2-1 


E fr 


(27) 


The discrete Fourier transform can be calculated over the frequency scale f t 
1)/A by the formula 


JV-l 

Ui = E u je 2i7rfit *A 

j=o 


and the discrete power spectrum is defined as 

/.' = //-. 


= * - 


(28) 


(29) 


Under these definitions, we have 


N 


/ 2 


U 


= y Ef /Delta 


(30) 


and the integral scale can be calculated as 


7 \i ,i n i — lim 


E 


o 4 ii r2 


(31) 


(see Hinze [7], p.65). 


Appendix 2 — Dominant Scales 


An alternative characteristic time scale based on the Fourier spectrum is the domi- 
nant (most energetic) spectral range. The spectral energy content in a logarithmic 
frequency band around / is estimated by the integral 


A E F (f)= [ [ E F (f) df = [ f E F {f)d{\og f)i 

J JAf J A log / 


(32) 


Accordingly, dominant scales are identified as the maxima of f Ep(f). Let us denote 
the dominant frequency as 
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When two signals are available the cross spectrum can be computed. The spectral 
distribution of u-v-correlations is obtained from the corresponding Parseval formula 
(of which Eq.(24) is a particular case): 


/ OO POO POO 

u(t)v*(t)dt = - ii(f)v*(f) f = ~ f ii(f)v*(f) dlog(f) 

-oo J — oo J —oo 

/ OO 

f R F (f)dlog(f). (33) 

-CO 

The asterisk denotes complex conjugation. The expression -fu(f)v*(f) = / Rp(f) is 
therefore a measure of the spectral distribution of Reynolds stress. Any peaks of this 
function identify dominant scales, as used in the text. 


Appendix 3 — Comparison of Integral and Domi- 
nant Scales 


This section addresses the differences between the integral and dominant scales de- 
fined above. There is no reason why the two methods would give the same result, but 
it is important that results can be compared. We will establish a basis for comparison 
by calculating the corresponding time/frequency scales from a simple expression of 
the power spectral density. We note that the formulae use wavenumber rather than 
frequency (to give the correct slope of the inertial spectrum), and that the spatial 
Fourier transform underlying these manipulations rely on the wavefunctions exp(i.Kx) 
rather than exp(2i.Tvkx). This difference with the formulae used for the frequency 
spectrum are consistent with Hinze’s notations, and results minor difference of ap- 
pearance between the spatial and temporal relations. For example, the length scale 
Lmax corresponding to a dominant wavenumber n. max is given by 


L 


max ,v max 


= l 


(34) 


to be contrasted with Eq.(8). 

A variant on the model of Pao (Hinze [7], p 240) will be used here, modifying 
a 3-D Kolmogorov spectrum to account for large- and small-scale effects. The 3-D 
energy spectrum we will use is 

E 3 (k) = as 2 / 3 k~ 5 / 3 exp ( — (/w/) 4 / 3 ) exp ( — (kL)~ 4 ^ 3 ) (35) 

where L is the size of the large eddies, r/ is the Kolmogorov microscale, and a is the 
Kolmogorov constant. A representative 3-D spectrum calculated from this formula is 
shown on Fig. 20. 
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Figure 20: Model of 3-D power spectral density for Re = 300. 


The parameters s, u, v and r/ are related through this .definitions 



CO 

II 

CO 

(36) 

and 

CO 

II 

(37) 

Accordingly: 



Es(k.) 

= au' 2 {n.L) 2/,3 k 1 exp ( — (/w/) 4 / 3 ) exp ( — {kL) 4 / 3 ). 

(38) 


Based on the formulae for isotropic spectra (Hinze [7], p 209) we calculated the 
equivalent 1-D spectrum Ei, in which directional aliasing results in constant value of 
Ei at low frequencies: 


Ei(ki ) 


= au 


n 


= jr hhu - f 

JOO (nL)- 2 l z n ~ 1 exp-(Kij ) 4 / 3 exp-(KL )~ 4 / 3 ^ _ tflyi K 


(39) 


The resulting 1-D spectrum is shown on Fig. 21. The normalization condition 


u 


a 

i 



(40) 


yields a value of the Kolmogorov constant a. 

Two degrees of freedom are associated with these spectra: we chose to fix L and 
vary Re = u' L/v as the independent parameter. 
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El Pa c 

Figure 21: Model of 1-D po wer spectral density for Re = 300. 

The estimates of the integral length scale L' (see Hinze [7], Eq. 3-74): 

_ t^E 1 { 0 ) 

mt ~ 2u'\ 

and of the most-energetic scale L max dehned by the maximum of k.\E\{k.\) were 
obtained numerically. The ratio of these length scales is equal to the inverse ratio of 
frequencies used on the data, since a common convection velocity is used (Taylor’s 
hypothesis) to convert from length- to time-scale. For a range of turbulence Reynolds 
numbers, the integral and dominant scales and their ratio are shown on Table 4. 



Table 4: Comparison of integral and dominant length scales 


Re 

Lint/L 

Lmax f L 

Lint/ L max 

30 

0.99 

1.60 

.62 

100 

1.04 

1.50 

.69 

300 

1.06 

1.45 

.73 

1000 

1.06 

1.42 

.75 

3000 

1.07 

1.41 

.76 

10000 

1.07 

1.45 

.74 

30000 

1.07 

1.45 

.74 

100000 

1.07 

1.45 

.74 


’or a model spectrum. 


Thus, under Taylor’s hypothesis, the characteristic frequencies 
integral scales would be about 30% as large as frequencies of the most 


associated with 
energetic range, 
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and the ratio can be expected to be almost independent of the Reynolds number of 
the turbulence. Therefore, we conclude that for an isotropic turbulence the values 
of the two large scales are comparable; while referring to different physical content 
(definitions), they characterize the same dynamical range in a way that is weakly 
dependent on the Reynolds number of the, turbulence 3 . However, the question of 
which scale is preferable will arise when we analyze the experimental data below. 


Appendix 4 — Wavelet Transforms and Spectra 

A wavelet g(t) is a function that meets the admissibility criterion [2], The wavelet 
transform u g of a signal u(t) is 

/ CO 

K 1 ' 2 u(T)g*(K.(T — t ) )cIt. (42) 

-oo 

The localization of the wavelet is governed by the variable t, and its duration by « _1 . 
The significant contributions to the integral is from the regions near the center of the 
wavelet. The magnitude of the wavelet transform u g (K,t) reflects how the shape of 
the wavelet matches the local shape of the signal, other things being equal. 

In this study, we used two types of wavelets: 

1. The Morlet wavelet is defined by the relation 

</>(,?’, ~ 0 ) = (cos 2? nr + i sin27r.r). exp ( — 2,r 2 7r 2 /~Q ) — exp ( — ~q/ 2 — 2,r 2 7r 2 /~Q), 

(«) 

i.e. a complex oscillation of unit frequency with a Gaussian envelope of half- 
width l./~o- The selection of ~ 0 reflects a compromise between time- and 
frequency-localization. A value between 5 and 10 is usually adopted, and we 
have used ~ 0 = 5 for all calculations reported below. The Morlet wavelet is 
shown on Fig. 22. 

2. The Mexican-hat wavelet (Mexha.t below) is distinguished by a single bump, 
represented analytically by 

t/ 2 (~) = (~ 2 — 1) exp ( — ~ 2 /2), (44) 

and also shown on Fig. 22. It provides better temporal localization than the 
Morlet wavelet, at the expense of the frequency resolution. 


3 In the case of the GE data, the relevant factor in this comparison will be the ‘effective Reynolds 
number’, with the frequency cut-off of the data playing the role of the Kolmogorov microscale. For 
example, with a ratio of the order of 20 between the energetic eddies and the cut-off, the effective 
Reynolds number would be of the order of 20 4 / 3 or about. 50. 
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Dimensionless period 



) i me ils i o r I e s s d u re ti o n 


Figure 22: Top: Real (solid line) and imaginary (dashed line) parts of the Morlet 
wavelet for ~ 0 = 5. Bottom: The Mexican hat wavelet. 
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Just as in the case of the Fourier power spectrum, the temporal/spectral distri- 
bution of energy formulae rest on Parseval’s theorem [2], 

/ oo 2 1 f °° f °° 2 f °° 

| u(t) | dt = — / / dt | u g (K,t) | = / 2Emo,Mx{k) (Ik, (45) 

-oo Cg JO J — oo JO 

by which the expression | u 2 {K,t) | 2 /2c 3 can be interpreted as the local energy spec- 
tral density. Em 0 and i.\i. are the Morlet and Mexhat wavelet alternatives to Ep 
describing the spectral distribution of energy 4 . 

Spectral decomposition of the Reynolds stress can also be achieved with Mexhat 
wavelets: it was shown in [9] that a simple inverse transform formula holds: 

r oo 

u(t) = / U 2 {k, t) da, (46) 

Jo 

where 

U 2 (a,t) = 7r- 1/2 2- 7/4 «- 1/2 i< 4 («/ J2,t). (47) 

Here u 4 is the wavelet transform of u by the wavelet g±{t) = d 2 (g 2 (t)) / dt 2 . Now 
consider the u — ^correlation. Using Eq.(46), it can be rewritten in the form [10] 


— uv = — 


coo r oo 


' — oo J — oo 

r oo r oo 


U 2 {K. u ,t)V 2 {K. v ,t)K. u K. v d{ log K u )d{ log K v ) 

/■' u. /■' u /»' l/. ( /■' u. • /■' !'/i lou /■• „ !</( lou P !. 


' — oo «/ — OO 


(48) 


Therefore, the Reynolds stress can be decomposed into contributions from the plane 
of the logarithms of durations of each component signal. A diagonal section (k u = k v ) 
is used in the next sections. 


4 For the Mexhat wavelet, the relationship between the Fourier frequency / and the inverse wavelet 
duration k - 1 was established [8] as f / n = v 2 2 . 5 / 2 7r = .2516, or approximately 4 wavelet durations 
per period of the corresponding sine wave. This relation is used to plot the wavelet power spectra 
in all figures below in terms of conventional frequency instead of wavelet duration scales. 

Note that the Morlet wavelet is defined in this report with a unit frequency of the oscillation 
rather than of the envelope, so that no further conversion is required. 
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